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ABSTRACT 

Modern sonars are large in size. This is a result of several factors 
connected with operation at low frequencies. Non-linear interaction of 
two high-frequency sound beams in water may provide a source of low- 
frequency sound without the use of large transducers . The theory of 
interaction of two plane waves is reviewed, and conversion efficiency 
is studied. Early experimental work is also reviewed. To provide new 
experimental data, primary sources of 385 kHz and 435 kHz are mounted 
side-by-side in an anechoic tank. The sound field between the parallel 
beams is explored. A difference-frequency signal is detected, but its 
amplitude is found to be greater than the theoretical value. Non- 
linearities in the receiving hydrophone are suspected as contributing to 
the 50-kHz signal level. Various techniques for eliminating psuedo- 
sound are discussed, and recommendations for further study are 
presented. 
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I. INTRODUCTION 



The history of surface-ship sonars has been characterized by ever- 
increasing size. This is a direct result of the lowering of operating 
frequency to take advantage of reduced attentuation in the medium. As 
signal frequency is lowered and wave length is increased, the transducer 
elements become large in order to maintain dimensions comparable to the 
wavelength. The lowering of sonar frequencies has also aggravated the 
problem of cavitation limiting. Since the cavitation threshold is lowered 
as frequency is lowered, the power that can be transmitted through a 
given area is diminished. Figure 1 shows that for 20 kHz, this threshold 
is only about 2 W/cm 2 . An increase in pulse length also lowers the 
threshold as shown in Fig. 2. Since mechanical steering of large trans- 
ducers is impractical, electrical steering of large arrays to form beams 
has to be employed. This further aggravates the cavitation problem by 
producing large pressure peaks relative to the average distribution over 
the array. 

If attenuation were not a limiting factor, sonar operation at frequen- 
cies of a few hundred kilohertz would be attractive. Shorter wavelengths 
would permit use of smaller elements, and as seen in Fig. 1, cavitation 
thresholds might be raised by a factor of ten or more. The resulting 
smaller transducer might be mechanically steerable, eliminating the 
pressure peaking of electrical steering. A reduction in transducer size 
would produce obvious weight and structural advantages in a surface 
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ship and would probably lower costs as well. In some applications, 
such as in active sonobuoys or acoustic torpedoes, small transducers 
are mandatory. An increase of operating frequency there could lead to 
smaller overall size, or transformation of space to power systems use, 
etc . 

The non-linear interaction of two high-frequency sound beams in 
water may be a useable source of low-frequency energy. The advantages 
of high-frequency operation as well as those of low frequency operation 
could be exploited simultaneously. In this process, stresses are 
established because of the different instantaneous particle velocities 
associated with the two waves. These stresses act to produce an 
acoustic source at the difference frequency. Hence, small transducers 
might be used to produce two high-frequency signals which would act to 
produce a low-frequency signal. The latter signal would propagate with- 
out suffering the large attenuation which limits high-frequency operation. 

Although a considerable amount of theoretical work has been done in 
this area, the experimental work which has been reported is limited in 
scope. Most of the work which has been done in water has been at 
frequencies above one megahertz. While beamwidth has been studied 
closely, little data exist on the absolute amplitude of the difference- 
frequency signals. The results of some work are clouded by psuedo- 
sound, which is developed due to non-linear summing at the hydrophone. 
The typical response of a non-linear device will contain functions of the 
second power. If the hydrophone is non-linear, the radiation pressures 
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of two incident waves may produce a signal with frequencies including 
the difference frequency. Separation of psuedo-sound and the true 
difference frequency signal has seldom been accomplished by direct 
experimental means . 

Most proposals for use of the low-frequency signal produced by high- 
frequency sources have concerned narrow-beamwidth sonar, frequency - 
scanning sonars, or wide-bandwidth applications. Parametric amplifiers 
have also been considered. References 1-4 suggest some of the interest- 
ing possibilities. 

The purpose of this work was to produce and measure a desired 
difference-frequency signal while using primary sources of a few hundred 
kilohertz. In pursuit of this goal, it was desired to determine the 
characteristics of psuedo-sound and to explore methods for its elimination 
in the problem . 
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II. THEORY 



The generation of sound by the non-linear interaction of two sound 
waves is best described by starting with the theory of sound generated 
aerodynamically . Lighthill, in Ref. 5, begins with the fundamental 
equations governing the propagation of sound. For a uniform non-viscous 
medium, without sources of matter or external forces, the linearized 
equations are: the equation of continuity. 



+ p -ini- = 

at Po ax 4 


o 

ii 



the equation of conservation of momentum. 



au t bp 

Do ‘ + * 

at ax 4 

and the equation of state. 


= o ; ( 2 ) 


P = D c 0 2 ; 


(3) 



where o / u, c 0 and p are the acoustic parameters. The resulting wave 
equation is 

&g P - c 0 2 7 s o = 0 . (4) 

at 3 

If externally applied fluctuating stresses are present, then a stress 
term must be added to Eq. (2), 



auj + a P 
M 2t ax t 

Equation (4) then becomes 


aT ti 

= • -sf • (5) 


S 2 p 2 3 

y? - - c “ 7 p = 


illu- . (6) 

ax. ax. 
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For the case of a fluctuating fluid flow within a larger body of fluid at 
rest, Lighthill shows that the stress term may be expressed as 

= 0 u i u j + Pi j c ° P . (7) 

Here, pu t Uj is the momentum flux tensor (Reynolds stress tensor) and 
p tJ is the compressive stress tensor. The last term represents the 
stresses of a simple hydrostatic pressure field. Viscous stresses and 
heat conduction are neglected. The fluctuating Reynolds stresses, 
corresponding to the variation in momentum flux across the fluid bound- 
aries, and the fluctuating applied pressures serve as sources of sound. 
Kinetic energy is converted to acoustic energy. Lighthill shows that T tJ 
is the strength per unit volume of a relatively weak quadrupole field. 

Westervelt, in Refs. 6 and 7, applies Lighthill' s work to the subject 
of two sound waves of different frequency, where Uj and Uj are taken as 
the instantaneous particle velocities of the two waves. Using Eq. (7), 
Westervelt rewrites Eq . (6) as 

□ 2 p = <p u ‘ u d • (8) 

He then shows that pressure p^ in the difference -frequency wave can be 
expressed as 



□ P c = ’ 

s 2c 0 






Here, p^ is given as a function of the primary-wave variables which are 
subscripted i. The substitutions 



Po C o 



and 7 2 Pj 3 = D 2 p, 2 + (p t s ) (10) 
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lead, finally, to the useful expression 



- Po 



aq 

bt 



( 11 ) 



where the source density function is 



is.* t 1 + 1 Cc 1 ir (p ‘ 2) • 



(12) 



-dp y p = 

For a given temperature, the expression in brackets is a constant. In 
water, its approximate value is three [Ref. 2], For convenience, it is 
taken as unity in what follows. 

If the primary wave is now taken as the sum of two colinear plane 
waves travelling to the right, acoustic pressure at a distance x may be 
taken as 

P 4 = Pi + P 2 = P i e Gl X cos (o\ t -k x x) + P 2 e a ^cos (o) 2 t-k^x) . 

(13) 



The variables are the usual ones. Squaring, 

p t 2 = ? x 2 e~ 2 a ' x C os» (aitr-kjx) + P 2 2 e -sa2X cos 2 (co^-k^x) + 

P,Po e [ cos [ ( Wl + co 2 ) t - (k x + k 2 )x] + 

cos [ («! - <o 2 )t - (k x - k 2 )x] } . (14) 

The first two terms lead to frequencies which are twice the primary 
frequencies. The last term contains both sum and difference frequencies. 
Following the development by Berktay [Ref. 2], the difference-frequency 
term is expressed as 

P 1 P 2 e ^ 0il + a 2 ^ X cos ( fit - K x) , (15) 

where fi = - w 2 and K = k 1 - k 2 . Then, from Eq. (12), the source 

density function is 
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q = P 4 P 1 P 2 e tol+(,z)x sin ( fit - Kx) . (16) 

Do C 0 

An elemental volume will contain a source of strength q(x,t) 6x6y6z, 
and the pressure at a distance r from the secondary source will be 



6 p (r,t) = - e ar 6x6y6z [q(t-“-) } , (17) 

s 4 77 r St c 0 

where a is the attenuation coefficient at the difference frequency. The 

total field at a point can be obtained by integrating the pressure field 

due to all such sources . 

Berktay treats the sources as filling thin wafers along the axis of 
the two primary beams . Each wafer contains sources of identical 
amplitude and phase. It is assumed that the difference-frequency waves 
propagate linearly and that they add in phase to the right. 

For the geometry of Fig. 3 and for R >>x, Berktay integrates Eq. (17) 
and finds the resulting pressure at the difference frequency: 

P s (R ' e) e ~ aR — ° - lpt - T - KR ~ 13 ) . (18) 

S 4 77 p 0 C 0 K r «2 . ro >7 .. .2 ® 



C A 3 + [2Ksin s f ] 2 H 



Symbols not heretofore defined are: 



S = cross-sectional area of the wafer; 
A~ a i + a 2 - a ; 

0 = arctan A/2K sin 2 ( -| ) ; 

0 = angle from the axis; 

R = distance from primary sources. 



Limiting analysis to the axis of the system, the pressure amplitude is 



P (R,0) 
s 



P.Po Q 2 S -aR 
477 Do c 0 4 RA 



(19) 
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Another convenient form, which follows immediately, is 



p (R o, -aR 

P s (R ' 0) - 2.C/RA 6 



( 20 ) 



Here, W x and W 2 are acoustic powers in the two primary waves. 

Summarizing, if two acoustic plane waves of different frequency 
propagate simultaneously along coincident paths, new waves are gener- 
ated due to non-linear interaction. The frequency of one of the new 
waves is equal to the difference frequency of the primary waves. Its 
pressure amplitude is proportional to the product of primary-^wave pres- 
sures and to the square of the difference frequency. Because of the 
quadrupole nature of the source, the difference-frequency signal is 
inherently weak. 

It should be noted that a similar treatment has been given to 
cylindrically spreading waves. The results are comparable and will not 
be covered here. The plane-wave case provides adequate background for 
basic experimental work. 
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